Expansion of a coherent array of Bose-Einstein condensates 
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We investigate the properties of a coherent array contain- 
ing about 200 Bose-Einstein condensates produced in a far 
detuned ID optical lattice. The density profile of the gas, im- 
aged after releasing the trap, provides information about the 
coherence of the ground-state wavefunction. The measured 
atomic distribution is characterized by interference peaks. 
The time evolution of the peaks, their relative population 
as well as the radial size of the expanding cloud are in good 
agreement with the predictions of theory. The 2D nature of 
the trapped condensates and the conditions required to ob- 
serve the effects of coherence are also discussed. 

PACS numbers: 03.75.Fi, 32.80. Pj 



Coherence is one of the most challenging features ex- 
hibited by Bose-Einstein condensates. On the one hand 
it underhes the superfluid phenomena exhibited by these 
cold atomic gases. On the other hand it characterizes in 
a unique way their matter wave nature at a macroscopic 
level. Coherence requires that the system be character- 
ized by a well defined phase, giving rise to interference 
phenomena. 

After the first interference measurements carried out 
on two expanding condensates at MIT the experimen- 
tal study of interference in Bose-Einstein condensed gases 
has become an important activity of research opening the 
new field of coherent atom-optics. 

The possibihty of confining Bose-Einstein condensates 
in optical lattices has opened further perspectives in the 
field 1^ . Bose-Einstein condensates confined in an optical 
standing-wave provide in fact a unique tool to test at 
a fundamental level the quantum properties of systems 
in periodic potentials. The observation of interference 
patterns produced by an array of condensates trapped 
in an optical lattice was already used as a probe of the 
phase properties of this system ||^,^ also allowing to proof 
the phase relation in an oscillating Josephson current . 
In the interference effect has been used to explore the 
emergence of number squeezed configurations in optically 
trapped condensates. 

The main purpose of this paper is to investigate the 
ground state properties of the system of a fully coherent 
array of condensates. To this aim we have explored the 
interference pattern in the expanded cloud, reflecting the 
initial geometry of the sample. 



The basic phenomenon we want to explore is the atom 
optical analog of light diffraction from a grating. The 
analogy is best understood considering a periodic and 
coherent array of identical condensates aligned along the 
X-axis. In momentum space the order parameter takes 
the form 

fc=0,±l..±fcM 

sin[(2/cM + l)pxd/2K\ 



saiPxd/2Ti 



(1) 



where k labels the different sites of the lattice, 2fcAf + 1 is 
the total number of sites (in the following we will assume 
kM ^1) and d is the distance between two consecutive 
condensates. The quantity nQ(px) = |^'o(pa;)P is the mo- 
mentum distribution of each condensate (see Eq. (||) be- 
low). The momentum distribution of the whole system, 
given by n{px) —\ '^{px) is affected in a profound way 
by the lattice structure and exhibits distinctive interfer- 
ence phenomena. Actually the effects of coherence are 
even more dramatic than in the case of two separated 
condensates Indeed, in the presence of the lattice the 
momentum distribution is characterized by sharp peaks 
at the values Px = n2TTh/d with n integer (positive or neg- 
ative) whose weight is modulated by the function no{px)- 
Furthermore, differently from the case of two separated 
condensates, interference fringes appear only if the ini- 
tial configuration is coherent. In principle the momen- 
tum distribution can be directly measured in situ using 
2-photon Bragg spectroscopy. This possibility has been 
already implemented experimentally for a single conden- 
sate 0. However, the very peculiar structure of (0) is 
expected to influence in a deep way also the expansion of 
the atomic cloud after the release of the trap. The width 
of the central peak (n = 0) of the momentum distribu- 
tion is of the order Apx ~ fi/Rx where Rx ~ kj^d is half 
of the length of the whole sample in the x-direction and 
the corresponding atomic motion, after the release of the 
trap, will be consequently slow. On the other hand the 
peaks with n ^ carry high momentum and the center 
of mass of these peaks will expand fast according to the 
asymptotic law 
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The occurrence of these peaks is the analog of multiple 
order interference fringes in light diffraction. 

We create an array of BECs of ^^Rb in the \F = 
l,mF = —1) state by superimposing the periodic opti- 
cal potential Vopt of a far detuned standing-wave on the 
harmonic potential Vb of the magnetic trap. For a more 
detailed description see |^,^. The resulting potential is 
given by 

V^Vb + Vopt 

= {ojIx"^ + ojliy^ + z^)) + sEr cos^{qx + |) (3) 

with m the atomic mass, ojx = 27r x 9 Hz and uj±^ — 
2tt X 92 Hz the axial and radial frequency of the mag- 
netic harmonic potential and x lying in the horizontal 
plane. In (^) s is a dimensionless factor, q — 2tt/\ is 
the wavevector of the laser light creating the standing 
wave and producing local minima in Vopt separated by 
d = A/2 and Er = n^q^/2m - 2TTh x 3.6 kHz is the re- 
coil energy of an atom absorbing one lattice photon. By 
varying the intensity of the laser beam (detuned 150 GHz 
to the blue of the Di transition at A = 795 nm) up to 
14 mW/mm^ we can vary the intensity factor s from 
to 5. We calibrated the optical potential measuring the 
Rabi frequency {VIr) of the Bragg transition between the 
momentum states —hq and +fiq induced by the standing 
wave. The intensity factor is then given by s = 2h^lR/ Er 

The procedure to load the condensate in the com- 
bined (magnetic -|-optical) trap is the following: we load 
^''Rb atoms in the magnetic trap and cool the sample via 
rf-forced evaporation until a significant fraction of con- 
densed atoms is produced. We then switch on the laser 
standing-wave and continue the evaporative cooling to a 
lower temperature (T^Tc). Typically, the BEG splits 
over ~ 200 wells, each containing 100 ~ 500 atoms. Af- 
ter switching off the combined potential we let the system 
expand and take an absorption image of the cloud at dif- 
ferent expansion times t^xp- 

In Fig. [T]A we show a typical image of the cloud taken 
at t 



exp 



29.5 ms, corresponding to a total number of 
atoms N — 20000 and to a laser intensity s — b. From 
the images taken after the expansion we can determine 
the relative population of the lateral peak with respect to 
the central one. The experimental results for the relative 
population of the first lateral peak as a function of the 
laser intensity s are shown in Fig. 2|. 

The structure of the observed density profiles is well 
reproduced by the free expansion of the ideal gas where 
the time evolution of the order parameter, in coordinate 
space, takes the form: 



*(x,t) = 



1 



(27r) 



ip^x/h —ip^t/2mh 



(4) 



For a realistic description of '^{px) we have improved 
the simple ansatz (^ in order to account for the k- 
dependence of the number of atoms Nk contained in 



each well. Due to magnetic trapping, the central con- 
densates with k <g; kM will be in fact more populated 
than the ones occupying the sites at the periphery. We 
have accounted for the modulation by the simple law 
Nk = No{l — k"^ /k\jY which will be derived below. For 
^'o we have made the Gaussian choice 



-^oiPx) oc exp[-p^cr72r] 



(5) 



corresponding, in coordinate space, to '^{j{x) ^ 
exp[— a;^/2cr^]. Using (j|) it is immediate to find that 
the relative population of the n ^ peaks with respect 
to the central one [n — 0) obeys the simple law 



Pn = exp[-47r^nV7d^] 



(6) 



holding also in the presence of a smooth modulation of 
the atomic occupation number Nk in each well. Result 
(^) shows that, if cr is much smaller than d the inten- 
sity of the lateral peaks will be high, with a consequent 
important layered structure in the density distribution 
of the expanding cloud. The value of cr, which charac- 
terizes the width of the condensates in each well, is de- 
termined, in first approximation, by the optical confine- 
ment. The simplest estimate is obtained by the harmonic 
expansion of the optical potential (^) around its minima: 
y = Efc(l/2)"^w^(a; - kdf with Cj^ = 2^sER/h, yield- 
ing (7 — d/ijTs^^^). However this estimate is not accurate 
except for very intense laser fields. A better value is 
obtained by numerical minimization of the energy using 
the potential (^ and the wavcfunction (||). This gives 
a/d — 0.30, 0.27, and 0.25 for s = 3, 4 and 5 respec- 
tively. The predicted results for the density distribution 
n{x) = |\E'(x)p evaluated for s = 5 and t = 29.5 ms are 
shown in Fig. |i|B (continuous line). 

From the above calculation we can also determine the 
relative population P„ of the n ~ \ peak as a function 
of the intensity factor s. This is shown in Fig. || to- 
gether with the experimental results. The good com- 
parison between experiment and theory reveals that the 
main features of the observed interference patterns are 
well described by this simple model. 

The ID model discussed above can be generalized to 
3-D through the ansatz 



*o(r) = 



*fc(r±) 



(7) 



fc=0.±l..±fcj 



which can be used, through a variational calculation, to 
describe the ground state of the system in the presence 
of the optical potential, magnetic trapping and two-body 
interactions. For sufficiently intense optical fields the 
value of a is not significantly affected by two-body in- 
teractions, nor by magnetic trapping. On the other hand 
interactions are important to fix the shape of the con- 
densate wave function in the radial direction. Neglecting 
the small overlap between condensates occupying differ- 
ent sites and using the Thomas-Fermi approximation to 
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determine the wave function in the radial direction we 
obtain the result 



9 



(8) 



where {R±)k — -\/2/ifc/maj^ is the radial size of the fc-th 
condensate, g depends on the scattering length a through 
the relation g — 'inh^a/m, while 



(9) 



plays the role of an effective fc-dependent chemical po- 
tential. The value of kM is fixed by the normalization 
condition N N^. and is given by 



- 

Km — 



2hLL! / 15 ^ a d 



2/5 



(10) 



In (|o|) u! — (ujxOJ^)^^^ is the geometrical average of the 
magnetic frequencies, aho — y^Ti/mZJ is the correspond- 
ing oscillator length and a is the s-wave scattering length. 
From the above equations one also obtains the result 
Nk = NM - k'^/kljf with A^o = (15/16)7V/fcM- Equa- 
tions (p[jlO|) generalize the well known Thomas-Fermi re- 
sults holding for magnetically trapped condensates |l^] 
to include the effects of the optical lattice. 

Neglecting two-body interaction terms in the determi- 
nation of the Gaussian width in the a;-direction is a good 
approximation only if /i^ is significantly smaller than the 
energy fiCox- This condition is rather well satisfied in the 
configurations of higher lattice potential employed in the 
experiment. For example, using the typical parameter 
iV = 5 X 10"' for the total number of atoms and the val- 
ues ZJ — 2tt X A2 Hz and a/uho — 3.2 x 10^'^, we find, 
for s = 4, (Da; ~ 27r X 14 kHz, ^k=o 27r?i x 0.5 kHz and 
/cm ~ 100, corresponding to Nq ~ 500. Notice that with 
these values the condition Hk S> ?ioj± required to apply 
the Thomas-Fermi approximation is rather well satisfied 
for the central wells. 

The fact that fj,k turns out to be significantly smaller 
than hibx not only explains why the interference patterns 
emerging during the expansion are well described by the 
ideal ID model for the array used above, but also points 
out the 2D nature of the condensates confined in each 
well. In this context it is worth pointing out that the 
bidimensionality of these condensates is ensured up to 
temperatures of the order oi kgT hoj^, which is sig- 
nificantly higher than the expected value of the critical 
temperature for Bose-Einstein condensation. Our sam- 
ple can then be used also to explore the consequence of 
the array geometry on the critical phenomena exhibited 
by these optically trapped Bose gases |12|. 

The above discussion permits also to explain the be- 
haviour of the radial expansion of the gas. In the pres- 
ence of the density oscillations produced by the optical 
lattice the problem is not trivial and should be solved 



numerically by integrating the GP equation. However, 
after the lateral peaks are formed, the density of the 
central peak expand smoothly according to the asymp- 
totic law R±{t) = R±{0)u!±texp, holding for a cigar 
configuration in the absence of the optical lattice p3[ . 
In Fig. ^ the linear law is plotted using the expression 
R±{0) '-^ {R±)k=o = kMduJx/uj± derivable from Eq. (||) 
for the condensate occupying the central well. This 
choice for -R_l(0) is justified if the population of the lat- 
eral interference peaks is small so that their creation does 
not affect the radial expansion of the system. 

Let us finally discuss the conditions required for 
our system to exhibit coherence. At zero tempera- 
ture the coherence between two consecutive conden- 
sates in the array is ensured if Ec «C Ej, where Ec 
and Ej are the parameters of the Josephson Hamil- 
tonian for two adjacent condensates jlj]. In particu- 
lar Ec — 2dfik/dNk is the interaction parameter while 

x=o is 

the Josephson parameter describing the tunneling rate 
through the barrier separating two consecutive wells. In 
our case (s = 4 and N = 5x 10*) we find Ec ^ 2Trh x 1 Hz 
for the most relevant central condensates {k <C fcjvf) 
while, by solving numerically the Schrodinger equation 
in the presence of the optical potential Kpt of Eq. (||), 
we find Ej ^ 2nh x 600 kHz. The value of Ej is so 
large that one can safely conclude that the ground state 
of the system is fully coherent and that the effects of the 
quantum fluctuations of the phase will be consequently 
negligible. This reflects the fact that, even for the largest 
values employed for the laser power the overlap between 
consecutive condensates is not small enough. The value 
of Ej is also much higher than the values of fcsT used 
in the experiment, so that also the effects of the ther- 
mal fluctuations of the phase of the condensate can be 
ignored. This suggests that the fringes associated with 
the expansion of the condensate will remain visible up to 
the highest values of T, corresponding to the critical tem- 
perature for Bose-Einstein condensation. We have car- 
ried out experiments at different values of T where the 
signal obtained by imaging the expanding cloud can be 
naturally decomposed in two parts: an incoherent com- 
ponent due to the thermal cloud which is parametrized 
by a classical Gaussian Boltzmann distribution, and a 
Bose-Einstein component exhibiting the interference ef- 
fects discussed above. In our experiment the interference 
peaks are visible up to fcsT ~ 27rh x 4.2 kHz. In order 
to point out the effects of the fluctuations of the phase 
one should lower the value of Ej by orders of magnitude. 
This can be achieved by increasing significantly the laser 
power generating the optical lattice. Such effects have 
been recently observed in the experiment of |^ . 

In conclusion, we have investigated the consequences of 
coherence on the properties of an array of Bose-Einstein 
condensates. We have observed peculiar interference pat- 
terns in the density of the expanded cloud, reflecting the 
new geometry of the sample and discussed on a theoret- 
ical basis some key features exhibited by these optically 
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trapped gases. Further studies in this direction include 
the possible effects of thermal decoherence ||l^,|l6| in the 
presence of tighter optical traps and the emergence of 2D 
effects in the thermodynamic properties of these novel 
systems. 
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FIG. 1. A) Absorption image of the density distribution of 
the expanded array of condensates. B) Experimental density 
profile (crosses) obtained from the absorption image (A) inte- 
grated along the vertical direction. The wings of the central 
peak result from a small thermal component. The continuous 
line corresponds to the calculated density profile for the ex- 
panded array of condensates for the experimental parameters 
(s = 5 and t^xp ~ 29.5 ms). 
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FIG. 2. Experimental (circles) and theoretical (triangles) 
values of the relative population of the n = 1 peak with re- 
spect to the n = central one as a function of the intensity 
factor s of the optical potential Vopt- 
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FIG. 3. Radial size of the central peak as a function of the 
expansion time. Experimental data point are compared with 
the expected asymptotic law J?x ~ R±{0)uj±texp- 
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